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ABSTRACT 

Algol (/? Per) is an extensively studied hierarchical triple system whose inner pair is a prototype 
semi-detached binary with mass transfer occurring from the sub-giant secondary to the main-sequence 
primary. We present here the results of our Algol observations made between 2006 and 2010 at the 
CHARA interferometer with the Michigan Infrared Combiner in the H band. The use of four telescopes 
with long baselines allows us to achieve better than 0.5 mas resolution and to unambiguously resolve 
the three stars. The inner and outer orbital elements, as well as the angular sizes and mass ratios for 
the three components are determined independently from previous studies. We report a significantly 
improved orbit for the inner stellar pair with the consequence of a 15% change in the primary mass 
compared to previous studies. We also determine the mutual inclination of the orbits to be much closer 
to perpendicularity than previously established. State-of-the-art image reconstruction algorithms are 
used to image the full triple system. In particular an image sequence of 55 distinct phases of the 
inner pair orbit is reconstructed, clearly showing the Roche-lobe-filling secondary revolving around 
the primary, with several epochs corresponding to the primary and secondary eclipses. 
Subject headings: binaries:eclipsing - infraredistars - stars: imaging - stars: individual (Algol, Bet Per, 
HD 19356) - techniques: interferometric - techniques: image processing - Online- 
only material: animation, color figures 



1. INTRODUCTION 

Algol (P Per, HD 19356) is one of the most extensively 
studied stellar systems in the history of astronomy. It 
is a triple hierarchical system where an inner eclipsing 
binary is orbited by a more distant outer star. The inner 
system of orbital period 2.87 days is a prototypical semi- 
detached binary in nearly circular orbit. The primary, 
Algol A, is a main-sequence star of type B8V. The sec- 
ondary, Algol B, is a subgiant of type K2IV. As Algol A 
and Algol B are thought to have formed roughly at the 
same time, one would expect the more massive star of the 
pair to become a subgiant first. However this supposi- 
tion is contradicted by all observations, showing Algol B 
to be lighter than Algol A. The solution to this so called 
"Algol paradox" is that most of the original mass from 
Algol B, which fills its Roche lobe, is thought to have 
been accreted by Algol A. As the inner pair is tidally 
locked, this creates a strong magnetic dynamo effect and 
non-thermal emission ob served in X-ray and r adio. Us- 
ing radio interferometry, iLestrade et al.l ([1993) detected 
positional displacement during the orbital revolution of 
the inner pair and identified Algol B as the source of 
radio emission. In X-rays, giant flares have been de- 
tecte d on Algol B (IWhite et al.l 119861 : iSchmitt fc Favatal 
1998: lSchmitt et al. 2003). Their sizes were estimated to 
a few tenths of the Algol B radius, and they were shown 
to follow active and q uiescent periodic cycles of 49 days 
(Richard s et al.l [2003). Ha profiles provide evidence for 
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mass transfer through a stream of gas fro m the Roche 
lobe of Algol B to Algol A (Gillet et al. 1989 }). In parallel, 
observations with EX OS AT (Favata et al. 2000), XMM 
(jSchmitt fc Nessll2004f ) and Chandra ( Drake! l2QQ3h have 
shown evidence of an extended corona around Algol B. 
The presence of a large, permanent coronal loop around 
Algol B, oriented toward Algol A, has ver y recently been 
confi rmed by radio images from the HSA (jPeterson et al.l 
l2010h . 

The tertiary component, Algol C, wa s first discov- 
ered by radial velocity measurement (jCurtissI [1908) 
and revolves around the inner pair over a period of 
about 680 days without eclipsing. It is an Am star 
with metal and hydrogen lines presenting F1V char- 
acteristics (jRichardsl 1 19931 ). Earlier studies based on 
astrometric, photometric, and radial velocity data al- 
lowed the determinati o n of the masses o f the inner 
stars dHill et al.l 119711: iTomkin fc Lambertl Il978[ ) and 
of Algol C (jBachmann fc Hershevl Il975h Then light 
curve meta-analyses combining pre vious visible, in- 
frared and ultraviolet obser vations (|Soderhjelml 1980; 
iRichards fc M ochnackil [l988[ ) produced estimates of the 
orbital elements for the inner and outer orbits, as well 
as for the surface temperatures, logg, radii and masses 
of the three stars. With the advent of speckle interfer- 
ometry, it b ecame possible to spat i ally resolve the outer 
comp onent (jLabevrie et al.l Il974t iMcAlister fc Degioial 
1979), leading to th e determination of orbital parameters 
for the outer orbit (jBonneaulll979l ). Optical interferom- 
etry then improved upon these results with the Mark 



2 



III interferometer (|Pan et al.l |1993[ ). However, due to 
180° ambiguities in early spec kle observations, all stud- 
ies based on speckle re sults ([Bonneaul 11979 ; Soderl ^elml 
1980; iPan et all Il993f ) produced ascending node esti- 
mates for the inner o r outer orbits in disagreement 
with the polarimetric ([Rudy fc Kempl Il978f ) and radio 
(|Lestrade et al.l [1993) observations. Therefore, as for 
given inclinations, the angle between the inner and outer 
orbits solely depends on the difference of the ascending 
nodes, this ambiguity implied that the orbits were found 
either roughly perpendicular or coplanar. As coplanarity 
was ruled out by the shallow ecl ipse depth of photometric 
observations (Soderhjelm 1980), ad- hoc corrections had 
to be made to the a scending nodes of speck le-based stud- 
ies. More recently, iCsizmadia et al.l ((2009) resolved the 
inner pair with long baseline interferometry at the Cen- 
ter for High Angular Resolution Astronomy (CHARA), 
using the CLASSIC instrument working in the near in- 
frared K s band, in combination with Very Long Baseline 
Interferometer radio observations. They determined the 
orbit of the inner binary to be in prograde rotation, in 
disagreement with the retrograde movem e nt found pre- 
viousl y in radio by iLestrade et al.l (j 19931 ) . IZavala et al.l 
(|2010l ) recently solved these long-standing issues, by si- 
multaneously resolving all three stellar components in 
the optical with the Navy Optical Interferometer (NOI). 
Using angular referencing, the most precise inner and 
outer orbits to date were derived, with positional er- 
rors on the stellar components of the order of a milli- 
arcsecond. The outer orbit was shown to be prograde and 
the inner orbit retrograde. The latest Very Long Baseline 
Array ra dio measurements are i n agreement w i th thi s ori- 
entation (Pe terson et al.ll2011f ). I Zavala et al.l (|201Q[ ) also 
presented the first reconstructed interferometric image of 
all three stars, though the short baseline lengths used at 
NOI precluded full separation of the inner pair compo- 
nents. 

In this paper, we present the results of our Algol obser- 
vations with the Michigan InfraRed Combiner (MIRC) in 
the H band, making use of the long baselines available at 
CHARA. In Section [3l the relative positions of the three 
components are determined by model-fitting for the 55 
epochs of our data set. New estimates of the orbital pa- 
rameters are then derived in Section fusing these stellar 
positions. Finally, images of the inner orbit are recon- 
structed for each epoch in Section[5l where we discuss the 
potential evidence for the phenomena detected by radio 
and X-ray instruments. 

2. OBSERVATIONS AND DATA REDUCTION 

Our observations of Algol were conducted between 
2006 and 2010 at the Georgia State University CHARA 
Array interferometer using the MIRC instrument. The 
CHARA Array, located on Mount Wilson, CA, consists 
of six 1 m telescopes and is the largest optical /IR interfer - 
ometer array in the world (|ten Bru mmelaar et al. 2005). 
Its 15 baselines range from 34 m to 331 m, providing res- 
olutions up to ~ 0.5 mas at H band and ~ 0.7 mas at K 
band. Thus the data collected constitute the highest an- 
gular resolution data available to date on Algol. MIRC 
was used to combine four CHARA telescopes together 
for true interferometric imaging in H band, providing 
six visibilities, four closure phases and four triple ampli- 
tudes simultaneously in eight narrow spectral bands for 



each snapshot (|Monnier et al.|[2006h . 

Our observations total 20 nights, making the accumu- 
lated interferometric data on Algol one of the largest 
datasets to date at CHARA. In particular, the data cov- 
ers most phases of the inner pair eclipse. A complete ob- 
serving log is listed in Table [TJ including our calibrators. 
The spacing of Algol observations is not regular: Algol 
was sometimes observed only as a backup target, when 
atmospheric conditions did not allow for other fainter 
targets. The duration of most of these observations gen- 
erally did not exceed 1 h. Two array configurations were 
used, S1-E1-W1-W2 and S2-E2-W1-W2. The full com- 
plement of fringes was always obtained except on 2009 
August 18 when bad seeing prevented the acquisition of 
the El fringes. Both configurations are well-adapted to 
imaging, with roughly equal Fourier coverage in all di- 
rections. Figure [1] presents the typical Fourier coverage 
achieved within one of our final data sets. 

The standard observing procedures were followed, and 
the data were reduced usin g the standard MIRC pipeline 
described in iMonnier et al.l ([20071 ). After frame co- 
adding, background subtraction and a Fourier transform 
of the raw data, fringe amplitudes and phases are used to 
form squared- visibilities and triple products. On all data 
prior to 2010, the photometric calibrations are estimated 
using shutter matrix measurements and partial beam 
chopping. Starting in 2010 data, the calibration relies on 
the newly available (and far superior) photometric chan- 
nel techniq ue, allowing flux measurement in parallel with 
the fringes (|Che et al 1120101 ). During the next step, cali- 
brators are used to calibrate the drifts in overall system 
response. Table [2] presents the sizes of the calibrators, 
as estimated from our ongoing study of "good" interfer- 
ometric calibrators and based on three independent pho- 
tometry methods, adopting different color-magnitude re- 
lations to compute the stellar sizes. It has to be noted 
here that 37 And is currently suspected to be a bi- 
nary, though with a flux ratio between components large 
enough not to significantly affect the calibration. Af- 
ter calibration, the pipeline outputs the final calibrated 
power spectra and bispectra into file s compliant with the 
OIFITS standard (jPauls et al.ll2005f ). 

Based on the orbital parameters obtained by previous 
studies, we expect significant stellar movement on short 
timescales. The orbital period of the inner pair is about 
2.87 days, meaning that Algol B is moving relatively to 
Algol A by about 0.07 mas in an hour. The movement 
of Algol C also impacts the visibilities through high or- 
der lobes which oscillate rapidly with its separation from 
the inner pair (ranging from 15 to 70 mas), though this 
phenomenon is somewhat mitigated by its much longer 
orbital period (about 680 days). 

While model-fitting can possibly be made to account 
for these time dependencies, image reconstruction pack- 
ages currently cannot. At best, an image reconstructed 
using our data accumulated over an hour would be 
blurred, and at worse it would not converge at all to- 
ward a meaningful solution. Thus, the data have to be 
split into shorter sequences of no more than 10 minutes, 
finding a compromise between reasonable Fourier cover- 
age and temporal consistency. This splitting is applied 
to about 60% of the nights and results in a total of 55 
different files, each of which contains about 240 power 
spectra and 160 bispectra, sufficient for imaging without 
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blurring. Figures [2] and [3] illustrate how the data accu- 
mulated on 2009 August 12 have been split into chunks. 

3. MODEL FITTING 

The purpose of model-fitting here is to derive the posi- 
tion, shape and flux parameters of all three Algol compo- 
nents. These parameters will then be used to determine 
the orbital parameters of the inner orbit (noted A-B) 
and outer orbit (noted AB-C), as well as to generate re- 
alistic prior images for later image reconstruction. Our 
general procedure is to minimize the residuals between 
the interferometric observables derived from our current 
model and the corresponding interferometric data, using 
a Levenberg-Marquardt algorithm. Two possible paths 
were initially examined: either directly fit the orbital 
parameters on the full data set, or estimate a separate 
set of parameters for every single epoch and then fit the 
orbital parameters using these positions. As Algol B 
is non-spherical (because it is thought to fill its Roche 
lobe), the two-dimensional projection of its sur face onto 
the o bserving plane varies from night to night (Ric hardsl 
1992). The "global fit" approach would require a com- 
plex three-dimensional model of the three stars, while 
the "night-to-night fit" approach we finally settled for is 
more flexible, as well as more manageable in terms of 
computing power and memory requirements. 

Each star is modeled here as an ellipse of uniform 
brightness distribution. As the angular sizes of the com- 
ponents are only few times the array resolution, the use of 
limb-darkening coefficients would have a nearly negligible 
impact on the visibilities. In fact, limb-darkened models 
are not strictly necessary for orbit fitting and even detri- 
mental as they increase the number of parameters to es- 
timate. Furthermore we expect image reconstruction to 
retrieve the true flux distribution. A major effect to take 
into account when modeling the system is the bandwidth 
smearing effect on Algol C, i.e. Algol C's contribution to 
the visibilities is affected by its large separation from the 
delay-tracking center (Algol A) and the width of each of 
the MIRC spectral channels. With this effect, the total 
number of parameters to estimate for each epoch is 23: 
the principal axes and angles for Algol B and Algol C 
relative to Algol A (the available interferometric data 
does not give access to absolute positions), the relative 
flux contributions of Algol A and Algol B to the total 
flux (normalized to unity), the shape parameters of the 
three elliptical disks (semimajor axis, ellipticity, and an- 
gular orientation), and the bandwidth smearing factors 
for Algol C (one for each spectral channel). 

The x 2 formed from power spectra and bispectra is 
known to be multimodal: as the Levenberg-Marquardt 
is a gradient algorithm, the parameters have to be ini- 
tialized close to their true value to prevent falling into 
local minima. Attempting to estimate the 23 unknowns 
simultaneously would only lead to inconsistent results. 
We determined that the positions of Algol B and Al- 
gol C relative to Algol A constitute the parameters with 
the most influence on the x 2 - Therefore our solution to 
this convergence problem involves simultaneous system- 
atic grid searches on both Algol B and Algol C Cartesian 
coordinates (relative to Algol A), coupled to Levenberg- 
Marquardt optimizations on the remaining parameters. 
In addition, some parameters can safely be fixed, based 
on the literature of the expected physical properties of 



the sys tem. As Algol A and B are thought to be tidally 
locked (Richards 1992), the inner orbit is posed circular. 
The principal axis of Algol B is supposed to be along the 
direction A-B, as expected for a Roche- lobe filling sec- 
ondary. Algol A and Algol C are both assumed perfectly 
circular, because they are not rapid rotators and are not 
filling their Roche lobes. With these assumptions the y 2 
fit on the remaining number of parameters is then well 
constrained as each data set contains between 120 and 
240 power spectra as well as between 60 and 160 bispec- 
tra (triple amplitudes and closure phases). Out of the 
eclipse periods, lower and higher bounds are imposed on 
each parameter solely to enforce the physical meaning 
of parameters. During the eclipses, the aspect ratios of 
Algol B are kept small (< 1.05). 

The fitting is done in two steps - first a coarse grid 
search on Algol B and Algol C positions, then a finer 
one. The coarse search uses a 10 x 10 grid with spacing 
0.1 mas aimed at determining the approximate positions 
of the components. For each epoch, the center of the Al- 
gol C grid is initialized at its expected posi tion based on 
the o rbital elements from the NOI results (Za vala et al.l 
|201Q[ ), while the center of the Algol B grid search is sim- 
ply set to the origin (0,0) (i.e., Algol A position). For 
a given night, this setup of the coarse grid search covers 
all the probable positions of the components. Then the 
finer grid search refines the estimates, using a 10 x 10 grid 
with spacing 0.02 mas around the most likely positions 
found during the coarse search. To insure convergence 
during both the coarse and fine searches, the fluxes and 
angular diameters of the star are always optimized first, 
while the secondary ellipticity and the bandwidth smear- 
ing factors are kept fixed to their initial estimates; then 
the constraints on these are relaxed and all parameters 
are re-optimized. 

To compute the visibility contribution of Algol C, 
an analytic expression of the complex visibilities of a 
bandwidth-smeared ellipse is directly employed. For the 
visibility contribution of the inner pair, the presence of 
epochs where one component eclipses the other prevents 
the use of analytic expressions. Algol A and Algol B are 
therefore modeled as ellipses on an image array of 64 by 
64 pixels with a pixel size of 0.05 mas. To prevent prob- 
lems linked to the discretization of the ellipse sizes with 
this pixellation, the edges of the ellipses are smoothed 
by a sigmoid function, so that the x 2 varies continuously 
with the size parameters. The complex visibilities are 
then computed by applying a Discrete Fourier Transform 
to the image. Finally the model power spectra and bis- 
pectra are derived from the visibilities, and the reduced 
X 2 is computed. 

Table [3] reports these reduced x 2 values for every night. 
Typically the agreement between the model and the data 
is good, especially considering the unrealistic assumption 
of uniform brightnesses: the values lie mostly within the 
1.0-6.0 range, with higher x 2 values encountered during 
the eclipse epochs where this model breaks down. Fig- 
ure H] presents the residuals of a typical fit, which show 
no noticeable trend. 

The full output of our model-fitting software consists 
in reduced x 2 maps as a function of Algol B and Algol C 
positions, as presented on Figure [5] The most probable 
positions of Algol B and Algol C correspond to the peak 
value of these maps. For each map, the contour at la 
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(68% confidence level) is computed and fitted by an el- 
lipse. The positions of the peaks for Algol B and Algol C 
as well as the parameters of the error ellipses (major and 
minor axes, angles) are recorded in Table [3l 

The fits also produce estimates for the angular diame- 
ters of all three stars. The estimates are (j)A = 0.88 ±0.05 
mas, <pB = 1-12 ± 0.07 mas, and <pc — 0.56 ± 0.10 mas. 
For Algol B, this angular diameter corresponds to the 
major axis of the ellipse (the aspect ratio of the Algol B 
ellipse changes with the phase of the inner orbit, be- 
tween 1.04 and 1.22 for our data sets, which we attribute 
the projection of the Roche lobe). 

Concerning Algol C, we note that the difficulty of esti- 
mating the bandwidth smearing factors is reflected here 
as a larger error bar associated with its angular diam- 
eter. Our angular diameters are not direc tly compara- 
ble wi th the published values in Table 1 of I Zavala et al.l 
(2010|, as these suffered from a conversion error (i.e. the 
parallax values in milliarsceconds and parsecs have been 
confused). Nevertheless, a s they were bas ed on the lin- 
ear diameters est imated inlRichardsl (pj)93), we settle for 
comparison with Ric hardsl (j!993f ) only. We will assume 
the most recent value published for the orbital paral- 
l ax, 34.7 ±0.6 mas, independen tly determined by opt ical 
(jZavala et al.l 120101 ) and radio (jPeterson et al.l [201 If ) in- 
terferometry, and more prec ise than the latest H ipparcos 
estimate of 36.17±1.40 mas (|van Leeuw en 2009). Finally 
the linear stellar radii computed with this parallax are 
given in Table HJ our results are consistent wit h the ra- 
dius e stimates derived by light curve analysis in Richardsl 
(1993), though they admittedly have greater error bars 
due to the night-by-night nature of our fits and the uni- 
form brightness assumption. 

4. ORBITAL SOLUTION 

Using the estimates of the relative positions of the com- 
ponents established in Section [3j the orbital parameters 
for the inner and outer orbits are derived. 

The orbital elements for both orbits and the mass ra- 
tio mA/rriB are fitted simultaneously, again using the 
Levenberg-Marquardt algorithm. The pa rameters are 
i nitia lized around the published values in I Zavala et al.l 
(|2010f ). for which the corresponding reduced % 2 is 
about 13.0. After convergence our final best x 2 is 0.94. 
In addition to these orbital parameters, Kepler's third 
law applied to both the inner and outer orbits allows us 
to derive the mass ratio mc/(rriA + tub) independently 
from parallax assumptions. The parameters correspond- 
ing to the best fit are given in Tables HHHJ with error bars 
corresponding to a 68% confidence level. The probability 
distributions of the parameters were obtained with the 
bootstrapping-with-replacement technique (using 50000 
bootstrap samples) and thus takes into account corre- 
lated errors. In the case where distributions were found 
to be roughly symmetrical about the peak values, only a 
single error bar is given. Figures [6] and [71 show, respec- 
tively, the inner and outer orbits corresponding to our 
best fit, as well as the estimated positions of Algol B and 
Algol C as determined in the previous section. 

Most of our estimated orbital elements are in agree- 
ment _witlipreyjmis results based on light curve analy- 
sis ([R ichards 1993J) and radio and optica l inter ferometry 
(jZavala et al.l 120101 : iPeterson et al.ll2011f ). In particular 
the respective prograde and retrograde movements of the 



outer and inner orbits are unambiguously verified (see 
also the reconstructed image sequence of the inner orbit 
in the next section). The mutual inclination angle i m 
between both orbits is given by: 

cos im = cos i\ cos %2 ± sin i\ sin %2 cos(^i — ^2) (1) 

where i\ and %2 are the inclination of the inner and 
outer orbits, and ^1 and ^2 are the associated as- 
cending node angles. The mutual inclination deter- 
mines the period of orbital precession between the or- 
bits, and consequently the amplitude of the variations of 
i nclination of the close binary due to dynamica l effects 
(jBorkovits et al.(2004) Borkovitsl iForgacs-DaikaT). Pre- 
vious ly published values close to 100° (jKiseleva et al.l 
1998) would imply a detectable change of inclination 
of the close binary from 1° to 3° in the last century, 
which in turn would lead to observable variations in the 
minimum depth of the light curves of the eclipse and 
the disappearance of t hese eclipses after a few centuries 
(jCsizmadia et al.l l2009h . This is however contrary to 
the actual observations: Algol eclipses have been known 
since Antiquity and therefore the minimum depth is not 
thought to have significantly changed within error bars. 
The mutual inclination estimate derived from our boot- 
strap results is i m = 90°. 20 ±0°. 32. It is completely con- 
sistent with photometric observations, and much closer 
to exact per pendicularity than the most recent estimate 
of 95° ± 5° (jZavala et al.ll20Toh . 

The main point of disagreement between previous 
studies and our work is that we found the semi-major 
axis of the inner orbit to be slightly smaller, in turn im- 
plying a smaller value for the total mass of the inner 
binary. However, to disregard this shorter axis result we 
would have to invoke a 15% error in our estimation of 
the position of the center of mass of the secondary. As 
our analysis is essentially based on the determination of 
photocenters, a potential source of error to consider is 
the asymmetric brightness distribution on the surface of 
a Roche- lobe-filling star. We quantified the magnitude 
of such an effect, by simulating the equatorial brightness 
profile of Algol B, taking into account both gravity dark- 
ening and limb darkening as shown in Figure [H The 
positions of the photocenter and of the center of mass 
were found to coincide to 1%. Arguably the existence of 
an additional strong proximity effect (due to the primary 
heating of the atmosphere of the secondary) could also 
explain a photocenter shift. Detailed modeling of such 
effect however goes beyond the scope of this paper, re- 
quiring complementary data (imaging in Section [5] based 
on current data did not allow detection of the proximity 
effect). Finally, the semimajor axis of the inner binary 
was also found to be 15% smaller than prev ious optical 
results by rec ent radio (jPeterson et aT1l2011l ) and X-ray 
observations (jChung et al.ll2004[ ) , a result strikingly simi- 
lar to ours. While both these papers hypothesized the ex- 
istence of a coronal component around Algol B, no coro- 
nal signature has been detected in the near-infrared yet. 
Our result thus suggests the possibility that the center of 
light and center of mass actually roughly coincide, with 
previous optical results suffering from imprecisions. Our 
observations indeed benefited from much higher angu- 
lar resolution than previous interferometric studies and 
from the rece nt disambiguation of the orbital angles from 
IZavala et all (j2010f ). 
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Consequently we recommend that our orbits (like 
all previously published orbits based on interferometric 
data) should be interpreted as photocenter orbits. We 
also note that we did not notice significant differences 
between the center of lights estimated by model-fitting 
uniform ellipses and estimated by using the real bright- 
ness distribution. In Section [5] we demonstrate that the 
semimajor axes estimated from the disk model and from 
the reconstructed images differ by no more than 5%. 

The inner binary mass ratio m^/m^ = 4.56 ± 0.34 
agrees with the most recent radial velocity measur ements 
(jHill et al.l[l993l fl97ll : iTomkin fc Lambertl fl978h which 
measured tra/tub = 4.6±0.1. Our larger error bar is due 
to the lesser intrinsic precision of interferometry in this 
parameter, as the mass ratio influences the position of 
Algol C only by a small amount. Figure [9] shows that the 
X 2 of the Algol C orbit fit possesses a very flat minimum 
valley, resulting in the large error on the mass ratio. 

With the assumption of a parallax of 34.7 ± 0.6 mas 
(j Zavala et al.l 120101 : iPeterson et al.l l2011f ) , the masses of 
the three stars are completely determined and we re- 
port the results in T able [H Algo l A an d Al gol B are 
lighter than found by IZavala et al.l (j2010f ) and iRichardsl 
(199 3]). We note he re that the tertiary mass published in 
I Zavala et al.l (j2010h is numerically inconsistent with the 
corresponding semi-axes and periods, which seems to be 
due to slightly different parallax values adopted for the 
inner and outer orbit fits (C. Hummel, private communi- 
cation, 2011). With our hierarchical fit, we find the mass 
of Alg ol C to be larger than that published in Ric hards! 
(1993). Overall our mass estimates are consistent with 
all measurements (radial velocity and system dynamics) 
and give a sensibly different view of the system. In par- 
ticular the picture of the mass transfer within the inner 
binary may have to be modified. 

5. IMAGE RECONSTRUCTION 

The amount of data present in each MIRC data sets 
corresponds to a dozen minutes of observation each night. 
All data sets except one (2009 August 18) contain four- 
telescope data. Consequently there is enough phase 
information to attempt "model-independent image re- 
construction" on each night, where the prefix "model- 
independent" underlines that the method does not rely 
on a specific astrophysical model. 

The image reconstruction procedure belongs to the 
class of "ill-posed" inverse problems, as it attempts to 
reconstruct an image conventionally consisting in a large 
number of pixels (typically a few thousand) using the 
less numerous interferometric data points (typically a few 
hundred). Expressing the problem in a Bayesian frame- 
work shows that the solution can be given by the regu- 
larized maximum likelihood m ethod (jBaron et al.l 1201(1 
Thiebaut & Giovannelli 2010). The target image is the 
array of fluxes x = {xo, . . . # n -i} that minimizes the sum 
of a term linked to the data (the likelihood) and a term 
reflecting all prior information (the regularization) under 
the two constraints of image positivity and of normaliza- 
tion of the image to unity, i.e. in the mathematical form: 

x = argmin [x 2 ( x ) + ^R(x)} . (2) 

subject to positivity (Vi, Xi > 0) and normalization 
(X^r^Q 1 x i ~ -0- The likelihood term measures the dis- 



tance between the observed power spectrum and bispec- 
trum and the same quantities derived from the tentative 
image. Minimizing this term enforces the presence of flux 
as indicated by the data, but the minimization presents 
many local minima. The addition of the regularization 
function R(x) discriminates between these local minima, 
as well as prevents over-fitting, detrimental to image 
quality. The factor \i in Equation (j2j) controls the rela- 
tive weight of the y 2 and regularization terms, and can be 
chosen so that the actual reduced x 2 is roughly unity for 
the reconstructed image (the selection of the truly opti- 
mal ji is a difficult problem which goes beyond the scope 
of this paper). Most regularizers are computed in the 
image plane, but also possesses good frequency extrap- 
olation properties in the Fourier plane. Hence they al- 
low reconstructions to routinely achieve super-resolution, 
i.e., to produce images with an effective resolution typ- 
ically about four times greater than the physical array 
resolution. 

The imag e reconstruction software SQUEEZE 

(|Baron et al l |2010[ ) has been used to obtain the re- 
construction presented in this paper. It combines the 
Markov Chain Monte Carlo and gradient-descent ap- 
proaches used respectively by i ts predecessors MACIM 
([Ireland! 120061 ) and BSMEM jB aron fc Yound [2008). 
SQUEEZE offers a vast choice of regularizers, ranging 
from classic entropy to novel wavelet regularizations. 
As the images of Algol are expected to consist of 
three compact limb-darkened stars, we choose to use 
the total variation regularizer, an edge-preserving 
regularizer which ensures both the smoothness and the 
comp actness of the stars while strongly penalizing stray 
flux (jRudin et al.l 119921 ). In our implementation, the 
total variation is implemented as: 

R^x) = 5^-TV(x<) =^2\x i+1 -Xi\, (3) 

i i 

A second regularizer is also employed in the form of a 
prior image to initialize the Markov Chain to a sensible 
starting point as well as to prevent the exploration of 
pixels where the presence of flux is known to be very 
improbable. The prior image m is created by convolving 
the best model image of Section [3] to one fourth of the 
array resoluti on. The en tropic prior expression is the 
Burg entropy (Burg 1975): 

R2(x) = ~ logfa/rrii). (4) 

i 

The number of pixel elements in the Markov Chain is 
set to 4000, and the length of the Chain to 10000 itera- 
tions. A common size and pixellation is chosen for all 55 
images to allow easier comparison. The highest instru- 
mental resolution is given by the largest CHARA base- 
line (Sl-Wl, 330 m), corresponding to 0.5 mas in the 
H band. Taking into account the expected amount of 
super-resolution, our choice of pixellation for the recon- 
structed images is thus 0.1 mas. Meanwhile the distance 
from Algol C to the inner pair is at most 100 mas, thus 
requiring an image size of 1000 x 1000 pixels. This is 
a factor 100 more pixels than typical binary reconstruc- 
tions from non-simulated data in optical/infrared inter- 
ferometry (|Zhao et al.l 120081 ). Nevertheless all 55 recon- 
structions converged without notable issues, with final 
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reduced \ 2 between 1.0 and 2.0. 

Figure 1X771 presents the reconstruction for 2009 August 
12. The three stars are present within a small field of 16 
mas, Algol C being close to its periapsis, and they are 
well resolved. 

Algol A appears as a bright and near-perfect circular 
disk with a small amount of limb-darkening. Algol B is 
fainter and elongated toward Algol A, as expected due 
to it filling its Roche lobe. Algol C also appears elon- 
gated in the direction of Algol A-B, though this is only 
here a purely optical-numerical artifact due to bandwidth 
smearing. 

All the reconstructed images of the inner pair are pre- 
sented in Figure [TTJ with epochs sorted by increasing 
phase of the inner eclipse. A movie generated from this 
sequence is also available with the online version of this 
paper. These images should be interpreted with caution, 
due to the comparatively poor resolution of the inter- 
ferometer compared to the size of the stars, and to the 
typically limited dynamic range of image reconstruction 
with four telescopes (about 10:1). Overall the aspect of 
Algol A barely changes from image to image, while Al- 
gol B's elongation varies as expected with the phase of 
the eclipse. 

Due to the sporadic nature of our Algol observations, 
the phase coverage is unequal and incomplete, and in par- 
ticular there are no data available beyond phase 0.861. 
About 20% of the images correspond to snapshots taken 
during the eclipse epochs. Images with a phase below 
0.08 correspond to the primary eclipse, when Algol A 
is occulted by Algol B. The progressive separation of 
the stars is clear as the phase increases up to 0.25 (even 
though the data for these images were taken several years 
apart). Images with a phase between 0.45 and 0.55 cor- 
respond to the secondary eclipse. During part of the 
secondary eclipse, Algol B seems to completely disap- 
pears behind Algol A despite its larger size, an effect we 
attribute to its lower brightness and the low dynamic 
range. 

In a significant number of images the secondary is 
clearly limb brightened, a phenomenon we attribute to 
gravity darkening. We can neither rule out nor confirm 
the presence of spots on the s t ellar s urfaces, such as hy- 
pothesized in Richards (1990, Il992f ). Similarly, we are 
unable to confirm the existence of a stream of matter 
between the stars or of an active corona around Algol B 
such as detected in radio (while the 2009 August 18 im- 
age at phase 0.463 would seem to show some interesting 
interaction between the stellar pair, it is the least reliable 
of all images due to data availability on three telescopes 
only). 

Using our reconstructed images, we can now quantify 
the impact of our assumption of uniformly bright stars 
adopted in Section [3] As a function of the inner or- 
bit phase, Figure [12] presents the amount of correction 
to the principal axis of A-B introduced by using the 
brightness distributions given by the imaging. No sig- 
nificant trend is detected. The maximum correction is 
0.057 mas, roughly corresponding to half a pixel, and 
thus always within one standard deviation of the posi- 
tions given in Table [3j Consequently, this validates our 
modeling in Section [3] and increase our confidence in the 
derived orbital solution. 



6. CONCLUSION 

From the analysis of the CHARA/MIRC infrared data 
on Algol we derived new orbital elements of the triple 
system with unprecedented precision. The respectively 
prograde and retrograde nature of the outer and inner 
orbits is unambiguously confirmed and the mutual in- 
clination of the orbital planes is found to be extremely 
close to perfect perpendicularity. Our results suggest 
that the semi-major axis of the inner orbit is shorter 
than previously reported by other methods. While this 
effect could in theory be due to a shift of the photocen- 
ter with respect to the center of mass, we did not find 
evidence to support this assumption. If our orbits reflect 
the center of mass orbit, then the masses of the primary 
and secondary are found to be lighter than suggested by 
previous works. Using model-independent image recon- 
struction techniques, we reconstructed a sequence of 55 
images of the inner binary with an effective resolution of 
0.2 mas. This "movie" of Algol A and Algol B covers 
most phases of the inner pair, including the primary and 
secondary eclipses. While the array resolution did not 
allow us to confirm the detection of potential flares or 
spots on the surface of the stars, the Roche-lobe shape 
of the secondary is clearly visible. 

To further improve the orbital solution presented in 
this paper, three-dimensional and time-dependent im- 
age reconstruction algorithms are currently being devel- 
oped at the University of Michigan. By modeling the 
secondary as a Roche lobe and by imaging directly on 
spheroids, we expect to considerably enhance the qual- 
ity of the reconstructions, enabling us to detect potential 
surface features on the stars. 
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Table 1 

Days of observation of Algol with CHARA/MIRC 

UT Date Telescope configuration Calibrators 

2006 Oct 9 S2-E2-W1-W2 Zet Per 

2006 Oct 11 S2-E2-W1-W2 Zet Per 

2006 Oct 12 S2-E2-W1-W2 Zet Per 

2007 Oct 4 S1-E1-W1-W2 37 And, Zet Per 

2007 Nov 23 S1-E1-W1-W2 Gam Tri 

2008 Aug 18 S1-E1-W1-W2 Zet Per 
2008 Aug 19 S1-E1-W1-W2 Zet Per 
2008 Aug 20 S1-E1-W1-W2 37 And 

2008 Aug 21 S1-E1-W1-W2 37 And, Zet Per 

2009 Aug 10 S1-E1-W1-W2 Gam Tri 
2009 Aug 11 S1-E1-W1-W2 Gam Tri 
2009 Aug 12 S1-E1-W1-W2 Gam Tri 
2009 Aug 13 S1-E1-W1-W2 Gam Tri 
2009 Aug 18 S1-W1-W2 Gam Tri 
2009 Aug 19 S1-E1-W1-W2 Gam Tri 
2009 Aug 20 S2-E2-W1-W2 Gam Tri, Zet Per 
2009 Aug 21 S2-E2-W1-W2 Gam Tri 

2009 Aug 24 S2-E2-W1-W2 37 And 

2010 Aug 6 S2-E1-W1-W2 Gam Tri, 10 Aur 
2010 Aug 8 S2-E1-W1-W2 Gam Tri, 10 Aur 
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Figure 1. Typical uv coverage for one of our Algol "split" data set when using the telescope configuration S2-E2-W1-W2 (left) and 
S1-E1-W1-W2 (right). 



Table 2 

Calibrator sizes and la Error in Milli-arcseconds. 

Calibrator Uniform Disk Size (mas) Error (mas) Reference 

10 Aur 0.419 0.063 b 

37 And 0.682 0.030 a,b 

Gam Tri 0.522 0.033 a, b, c 

Zet Per 0.703 0.021 b 



References. — (a) iKervella et al.l (|2008l) ; (b) iBarnes et~aTl (|1978T) 
(c) iBonneau et all (|2006f) . 
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Figure 2. Squared visibilities collected on 2009 August 12, demonstrating how the dataset has been splitted into 5 epochs (2009 August 
12 a,b, . . . e). Both Algol visibilities and the uv coverage are significantly changing over the course of 20 minutes. For clarity we only show 
one-fifth of the full data set. 
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Figure 3. Closure phase data collected on 2009 August 12 presented as a function of wavelength and sorted by epochs. The error bars 
on the closure phases are of the order of 1°. 
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Figure 4. 

a data set 



Power spectrum residuals (left) and closure phase residuals (right) normalized by the data error, here for the 2009 August 12 
with a reduced x 2 of 1.16. 
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Algol B position - Likelihood map Algol C position - Likelihood map 
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Figure 5. Likelihood maps of Algol B (left) and Algol C (right) positions, for the typical 2009 August 12 a reconstruction. The solid line 
corresponds to the la confidence level, and the dashed line is the best-fitting error ellipse. 



Algol B position, relative to A at origin 

2l 




Figure 6. Orbit of the inner binary: (left) band of allowed orbits at 3cr for Algol B relative to Algol A (green), and best fit solution (solid 
black line) ; the error ellipses for each epoch are derived from the % 2 distribution maps obtained when fitting a uniform brightness model; 
(right) declination and right ascension as a function of the inner orbit phase, the zero phase corresponding to the light minima. 
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Algol C position, relative to the center of mass of A-B 
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Figure 7. Orbit of the outer binary: (left) band of allowed orbits at 3-cr for Algol C relative to the center of mass of the inner binary 
A+B (green), and best fit solution (solid black line); (right) declination and right ascension as a function of the outer orbit phase (the 
reference zero phase is at To given in Table [5}; the semi-major axes of the error ellipses are at worse 0.5 mas and thus barely visible. 




-2 2 4 

Equatorial position (R ) 

Figure 8. Simulated profile of the equatorial brightness of Algol B for a Roche lobe filling factor of 1.0. As the surface departs from 
sphericity, the radius-dependent gravity darkening effect induce a shift of the photocenter (dashed line) with respect to the center of mass 
(zero position). For this simulation, as we move along the equator, we compute the distance of the surface to the center of mass, using 
the Roche-von Zeipel equation. We then derive the corresponding local gravity and local temperature values. The brightness distribution 
is given as the combination of gravity darkening and limb darkening. We used the results from infrared light-curve fitting published in 
Richards (1990) to set the gravity darkening coefficient (f3 = 0.08, corresponding to the theoretical result for a convective envelope) and 
the quadratic limb-darkening coefficients (u\ = 0.680, U2 = —0.171). The reported brightness distribution is normalized to its maximum 
value. Overall we find that the position of the center of light differs from that of the center of mass by only 1%. 
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Algol triple reconstruction - 12 August 2009 
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Figure 10. Reconstructed image of the triple system for the 2009 August 12 a data set (final reduced x 2 =1.02). The three components 
are resolved: Algol A is nearly circular, Algol B is elongated as it fills its Roche lobe, and Algol C is elongated due to bandwidth smearing. 
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Figure 11. Evolution of the SQUEEZE reconstructions of the inner stellar system A-B with the phase of the primary eclipse. The phase 
is given as <E> = [(To — MJD)/Pab mod 1] with the orbital elements from Table [6] <E> = corresponds to the primary eclipse and the time 
of the light minima, and <I> = 0.5 corresponds to the secondary eclipse. 
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Figure 12. Difference between the principal axis of the inner pair estimated from the reconstructed images and estimated from uniform 
ellipse model-fitting, as a function of the inner orbit phase. The principal axis was defined as the vector joining the centers of light of 
Algol A and Algol B, Its computation was restricted to epochs where the stellar disks do not overlap. 
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Table 3 

Position estimates of Algol B and Algol C relative to Algol A 



Observation 


MJD 


X 2 /dof p B 


Ob 


Ca,B &b,B ^B PC Oc &a,C °b,C 


ip C 


(UT date) 


(days) 


(mas) 


(deg) 


(mas) (mas) (mas) (mas) (deg) (mas) (mas) 


(deg) 



2006 


Oct 9 a 


54017, 
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0, 


,103 


318.16 


44. 


,54 


-61. 


.91 


0, 


.217 


0, 


.140 


302.82 


2006 


Oct 9 b 


54017. 


.547162 


1 


,193 


1 


.565 


39.72 


0. 


,199 


0. 


,083 


310.04 


44. 


,55 


-61. 


,93 


0. 


.192 


0, 


,141 


299.08 


2006 


Oct 9 c 


54017. 


,558304 


2 


.778 


1 


.558 


38.24 


0. 


.170 


0. 


,094 


299.65 


44. 


,52 


-61. 


,88 


0, 


,229 


0. 


,158 


106.26 


2006 


Oct 11 a 


54019. 


.403382 


6 


.467 





.538 


109.70 


0. 


,286 


0. 


,129 


286.92 


45. 


,09 


-61. 


,73 


0. 


.331 


0, 


,202 


108.53 


2006 


Oct 11 b 


54019. 


.416411 


4 


,271 





.485 


87.79 


0. 


.283 


0. 


,238 


160.59 


45. 


,62 


-61. 


,73 


0. 


.421 


0. 


,179 


305.31 


2006 


Oct 11 c 


54019. 


,553074 


3 


,673 


1 


.033 


59.38 


0. 


.217 


0. 


,118 


227.40 


45. 


.74 


-61. 


.66 


0. 


,243 


0. 


.139 


119.75 


2006 


Oct 11 d 


54019. 


.555134 


3 


.188 





.983 


59.86 


0. 


.211 


0. 


,103 


222.20 


45. 


,88 


-61. 


.72 


0. 


.234 


0. 


.150 


298.48 


2006 


Oct 12 a 


54020. 


,441845 


1 


,124 


1 


.473 


35.46 


0. 


,271 


0. 


.183 


303.77 


46. 


,37 


-61. 


,28 


0. 


,309 


0. 


,184 


128.15 


2006 


Oct 12 b 


54020. 


,451949 


1 


.540 


1 


.396 


33.72 


0. 


.155 


0. 


,086 


307.40 


46. 


.35 


-61. 


,31 


0. 


.238 


0. 


,153 


136.16 


2007 


Oct 4 a 


54377. 


.358738 


2. 


,891 


1 


.914 


-137.73 


0. 


.555 


0. 


,134 


61.12 


18. 


,86 


108. 


,01 


0. 


,374 


0. 


.075 


67.59 


2007 


Oct 4 b 


54377. 


.372951 


2 


.098 


1 


.752 


-146.02 


0. 


,289 


0. 


.139 


64.24 


18. 


,81 


107. 


,87 


0. 


,330 


0. 


.100 


242.73 


2007 


Oct 4 c 


54377. 


,384093 


2 


,076 


1 


.670 


-146.21 


0. 


,367 


0. 


,207 


234.45 


18. 


.69 


108. 


,25 


0. 


,806 


0. 


.137 


64.15 


2007 


Nov 23 a 


54427. 


.234904 


3. 


,484 


2 


.200 


44.83 


0. 


,330 


0. 


.195 


103.13 


61. 


.06 


128. 


,69 


0. 


,508 


0. 


,157 


124.32 


2007 


Nov 23 b 


54427, 


,249729 


1. 


.975 


2 


.196 


43.80 


0. 


,252 


0. 


,159 


242.64 


60. 


.69 


127. 


,74 


0. 


,449 


0. 


,149 


125.02 


2007 


Nov 23 c 


54427, 


,265206 


2. 


.795 


2 


.198 


41.31 


0. 


,248 


0. 


,100 


232.43 


60. 


.74 


127, 


,75 


0. 


,254 


0. 


.150 


309.56 


2008 


Aug 18 a 


54696. 


,522101 


5. 


,493 


1 


.935 


48.69 


0. 


,234 


0. 


,139 


231.63 


44. 


,17 


-62, 


,15 


0. 


,344 


0. 


.109 


298.54 


2008 


Aug 18 b 


54696. 


,525448 


3. 


,876 


1 


.985 


49.72 


0. 


,326 


0. 


,158 


242.68 


43. 


.80 


-62 


,15 


0. 


,405 


0. 


,178 


301.15 


2008 


Aug 18 c 


54696. 


,536763 


4, 


,791 


1 


.984 


48.36 


0. 


,262 


0. 


,115 


226.02 


43. 


.99 


-62. 


,12 


0. 


.427 


0. 


.121 


293.89 


2008 


Aug 18 d 


54696. 


,539548 


3. 


,268 


1 


.984 


48.52 


0. 


,243 


0. 


,110 


231.84 


43. 


,83 


-62 


,15 


0. 


.401 


0. 


.121 


296.55 


2008 


Aug 19 


54697. 


,529045 


3. 


,965 





.534 


-92.89 


0. 


.219 


0. 


,130 


253.86 


44. 


.67 


-62 


,38 


0. 


.303 


0. 


.122 


116.32 


2008 


Aug 20 a 


54698. 


,455342 


4. 


,235 


1 


.829 


-144.44 


0. 


.318 


0. 


,214 


212.86 


44. 


,90 


-62 


,52 


0. 


,448 


0. 


,174 


128.07 


2008 


Aug 20 b 


54698. 


,463525 


1. 


,703 


1 


.730 


-146.33 


0. 


.278 


0. 


,264 


62.05 


45. 


.32 


-62 


,46 


0. 


,471 


0. 


,151 


127.81 


2008 


Aug 21 


54699. 


,475045 


3. 


,292 


2 


.065 


46.72 


0. 


,196 


0. 


,149 


273.80 


45. 


.59 


-61. 


,46 


0. 


,299 


0. 


,086 


305.89 


2009 


Aug 10 a 


55053. 


,531881 


4. 


,257 


1 


.908 


-131.04 


0. 


,276 


0. 


,175 


220.41 


15. 


,17 


102. 


,16 


0. 


,320 


0. 


,186 


131.00 


2009 


Aug 10 b 


55053. 


,503913 


3. 


,033 


1 


.818 


-130.08 


0. 


,302 


0. 


,153 


235.73 


15. 


.15 


102. 


,10 


0. 


.223 


0. 


,177 


318.71 


2009 


Aug 10 c 


55053. 


,525093 


4. 


,020 


1 


.877 


-131.65 


0. 


,258 


0. 


,179 


237.61 


15, 


.16 


102. 


,12 


0. 


,315 


0. 


,180 


122.35 


2009 


Aug 10 d 


55053. 


,528275 


4. 


,348 


1 


.884 


-131.65 


0. 


,279 


0. 


.182 


229.05 


15, 


,17 


102. 


,12 


0. 


.325 


0. 


,172 


130.78 


2009 


Aug 11 a 


55054. 


,497896 


3. 


,595 





.307 


107.63 


0. 


.261 


0. 


,167 


124.70 


16. 


.35 


102. 


,97 


0. 


,269 


0. 


,194 


118.47 


2009 


Aug 11 b 


55054. 


,513533 


7. 


.212 





.353 


94.87 


0. 


.194 


0. 


,153 


121.59 


16. 


.38 


102. 


,91 


0. 


,208 


0. 


.148 


310.20 


2009 


Aug 12 a 


55055. 


,480222 


1. 


,157 


1 


.728 


34.39 


0. 


.210 


0. 


,141 


211.54 


17. 


,18 


103. 


,62 


0. 


,124 


0. 


,053 


300.49 


2009 


Aug 12 b 


55055. 


,490959 


0. 


.724 


1 


.630 


34.83 


0. 


.231 


0. 


,163 


217.99 


17. 


,18 


103 


,68 


0. 


,191 


0. 


.155 


115.96 


2009 


Aug 12 c 


55055. 


,516855 


1. 


,169 


1 


.585 


33.70 


0. 


.211 


0. 


,100 


214.40 


17. 


.20 


103. 


,76 


0. 


,114 


0. 


,091 


289.18 


2009 


Aug 12 d 


55055. 


,497004 


0. 


,983 


1 


.639 


33.68 


0. 


,219 


0. 


.123 


222.33 


17. 


,19 


103. 


,69 


0. 


,099 


0. 


,093 


332.47 


2009 


Aug 12 e 


55055. 


,466371 


0. 


,578 


1 


.676 


35.63 


0. 


,194 


0. 


,071 


32.54 


17. 


.19 


103 


,65 


0. 


,115 


0. 


,062 


123.01 


2009 


Aug 13 a 


55056. 


,508177 


3. 


,446 


2 


.115 


-133.91 


0. 


,279 


0. 


,215 


267.54 


17. 


,79 


107. 


,45 


0. 


,266 


0. 


,188 


122.84 


2009 


Aug 13 b 


55056. 


,470177 


2. 


,552 


1 


.999 


-136.27 


0. 


,175 


0. 


,138 


-51.76 


17. 


.71 


107. 


,21 


0. 


,181 


0. 


,103 


299.71 


2009 


Aug 13 c 


55056. 


,463975 


1. 


,953 


2 


.027 


-136.49 


0. 


,207 


0. 


,137 


140.28 


17. 


.77 


107 


,42 


0. 


.153 


0. 


,099 


137.45 


2009 


Aug 18 


55061. 


,515703 


2. 


.238 





.794 


19.77 


0. 


,175 


0. 


,076 


143.17 


22. 


.44 


112. 


,29 


0. 


,360 


0. 


,062 


246.93 


2009 


Aug 19 a 


55062. 




5. 


A K 1 

,401 


2 


.111 


—136.85 


u. 


q a n 


u. 


,219 


lz4.Uo 


26. 


1 A 

. 14 


114. 


Q A 

,o4 


u. 


A QO 


u. 


.179 


126.58 


2009 


Aug 19 b 


55062. 


,497488 


2. 


.599 


2 


.083 


-138.78 


0. 


,369 


0. 


,200 


119.47 


23. 


,21 


114. 


.27 


0. 


,336 


0. 


,195 


129.47 


2009 


Aug 19 d 


55062. 


,503528 


3. 


,903 


2 


.076 


-139.32 


0. 


.229 


0. 


,137 


275.62 


23. 


,29 


114, 


,34 


0. 


.190 


0. 


,114 


302.48 


2009 


Aug 19 e 


55062. 


,506739 


7. 


,770 


2 


.020 


-141.33 


0. 


,087 


0. 


,068 


54.47 


23. 


,22 


114. 


.25 


0. 


.099 


0. 


.061 


122.80 


2009 


Aug 20 a 


55063. 


,416862 


2. 


,399 


1 


.516 


51.64 


0. 


,174 


0. 


,072 


141.62 


24, 


.38 


113, 


,59 


0. 


,195 


0. 


,129 


101.55 


2009 


Aug 20 b 


55063. 


,493563 


1. 


,571 


1 


.696 


49.13 


0. 


,178 


0. 


,090 


89.27 


24, 


,44 


113, 


,54 


0. 


,191 


0. 


,118 


239.39 


2009 


Aug 20 c 


55063. 


,503760 


1. 


,862 


1 


.731 


49.53 


0. 


,182 


0. 


,104 


90.96 


24. 


,44 


113, 


.48 


0. 


,180 


0. 


.162 


68.83 


2009 


Aug 21 a 


55064. 


,505563 


2. 


,316 





.206 


-47.38 


0. 


.240 


0. 


,134 


316.11 


25. 


,26 


115, 


,08 


0. 


,183 


0. 


,142 


106.54 


2009 


Aug 21 b 


55064. 


,516370 


4. 


,145 





.377 


-64.16 


0. 


.231 


0. 


,195 


330.11 


25. 


,26 


115, 


,07 


0. 


,171 


0. 


,131 


199.17 


2009 


Aug 24 b 


55067. 


,523787 


8. 


,119 





.717 


-109.20 


0. 


.367 


0. 


,175 


227.88 


28. 


.05 


117. 


.26 


0. 


.285 


0. 


.262 


291.42 


2010 


Aug 6 


55414, 


,467228 


1. 


.507 





.733 


-113.48 


0. 


.199 


0. 


,059 


65.96 


66. 


.78 


-56. 


.36 


0. 


,260 


0. 


,085 


133.57 


2010 


Aug 8 a 


55416. 


,489153 


0. 


.693 


2 


.135 


43.70 


0. 


,079 


0. 


,062 


403.26 


67. 


,73 


-55. 


,73 


0. 


.139 


0. 


.053 


128.93 


2010 


Aug 8 b 


55416. 


,495590 


0. 


,737 


2 


.130 


43.59 


0. 


,068 


0. 


,073 


64.08 


67. 


,73 


-55, 


,74 


0. 


,111 


0. 


,050 


126.83 


2010 


Aug 8 c 


55416. 


,499486 


0. 


.821 


2 


.144 


43.41 


0. 


,076 


0. 


,060 


120.88 


67. 


.69 


-55. 


,75 


0. 


,137 


0. 


.047 


125.75 


2010 


Aug 8 d 


55416. 


.501757 


0. 


,628 


2 


.149 


43.08 


0. 


,149 


0. 


,100 


312.13 


67. 


.72 


-55. 


,75 


0. 


,222 


0. 


,079 


132.75 



Note. — All angles are given east of north, ps and 6b are the polar coordinates of Algol B, while (T a ,B, <Jb,B-> and ipB are, respectively, 
the semimajor axis, semiminor axis, and angular inclination of its error ellipse. 



18 



Table 4 

Radius and Mass Estimates of Algol A, Algol B and Algol C. 





Richards et al. (1993) 


Zavala et al. (2010) 


This Work a 


Ra (Rq) 


2.90 ±0.04 




2.73 ±0.20 


Rb (Rq) 


3.5 ±0.1 




3.48 ±0.28 


Rc (Rq) 


1.7 




1.73 ±0.33 


M A (M ) 


3.7 ±0.3 


3.7±0.2 b 


3.17 ±0.21 


M B (M ) 


0.81 ±0.05 


0.8±0.1 b 


0.70 ±0.08 


M C (M ) 


1.6 ±0.1 


1.5±0.1 b 


1.76 ±0.15 



a As suming a parallax of 34.7 ± 0.6 mas (Zavala et al. 2010; Peterson et al. 
f20lTh . 

b Initially based on Richards (1993), then re-optimized using the Navy Optical 
Interferometer data. 



Table 5 

Orbital Parameters for the Outer Orbit. 



Orbital Element Zavala et al. (2010) This Work 



T (JD) 2446931.6 ±0.1 2446927.221^^ 

P (days) 679.85 ± 0.04 680.168 ± 0.046 

a (mas) 93.8 ± 0.2 93.43+g" H 

i (deg) 83.7 ± 0.1 83.66 ± 0.03 

Q (deg) 132.7 ± 0.1 132.66 ± 0.08 

e ' 0.225 ±0.005 0.227 ± 0.002 

uj (deg) 310.8 ± 0.1 310.02 ± 0.26 

m A /m B ... 4.56 ±0.34 

mC/(mA + mB) ... 0.456 ± 0.022 



Table 6 

Orbital Parameters for the Inner Orbit 



Orbital Element Zavala et al. (2010) This work 



T (JD) a 2441771.3395 2441771.353 ± 0.007 

P (days) 2.867328 2.867328 ± 5 x lO" 5 

a (mas) 2.3 ±0.1 2.15 ±0.05 

i (deg) 98.6 98.70 ±0.65 

Q (deg) 47.4 ± 5.2 43.43 ± 0.32 

e 
uj (deg) 



a Time of minimum light at the primary eclipse. 



